Next-to-leading QCD effect to the quark compositeness search at the LHC 
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We present the exact next-to-leading order (NLO) QCD corrections to the dijet production in- 
duced by the quark contact interactions at the CERN Large Hadron Collider (LHC). We show that 
as compared to the exact calculation, the scaled NLO QCD prediction adopted by the ATLAS Col- 
laboration has overestimated the new physics effect on some direct observables by more than 30% 
and renders a higher limit on the quark compositeness scale. The destructive contribution from the 
exact NLO correction will also lower the compositeness scale limit set by the CMS Collaboration. 
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The quark composite models have been studied exten- 
sively in the literature [H Q ■ It is assumed that quarks 
are composed of more fundamental particles with new 
strong interactions at a composite scale A, much greater 
than the quark masses. At energy well below A, quark 
contact interactions are induced by the underlying strong 
dynamics, and yield observable signals at hadron collid- 
ers. For example, the dijet production at the CERN 
Large Hadron Collider (LHC) could be largely modified. 
In the Standard Model (SM), the theory of Quantum 
Chromodynamics (QCD) predicts the jets in dijet events 
are preferably produced in large rapidity region, via small 
angle scatterings in t-channel processes. On the contrary, 
the dijet angular distribution induced by the quark con- 
tact interactions is expected to be more isotropic. The DO 
and CDF Collaborations at the Fermilab Tevatron have 
set limits on the scale A based on their dijet data Q. 
Recently, both the ATLAS and CMS collaborations have 
carried out similar analyses [3, 111 using the LHC dijet 
data with y/s = 7 TeV, in proton-proton collisions, and 
an integrated luminosity of about 3 pb^^. The observed 
limits for A at the 95% confidence level (CL) are 3.4 TeV 
and 4.0 TeV, respectively, which have already exceeded 
the previous limits. With more integrated luminosity col- 
lected, these limits will be further improved. 

The limits on the composite scale A were obtained by 
comparing the experimental dijet data with various the- 
ory predictions, including the SM next-to-lcading order 
(NLO) QCD corrections and contribution induced by the 
quark contact interactions which were handled differently 
in different experiments. While the CMS Collaboration 
included only the leading order (LO) contribution from 
the quark contact interactions Q, the DO, CDF and AT- 
LAS Collaborations included the "scaled NLO QCD cor- 
rection" which assumes the NLO correction (in terms of 
K-factors) to the dijet production from the contact inter- 
actions to be exactly the same as that from the SM QCD 
interactions [1, In this Letter, we present the exact 
NLO QCD correction to the dijet production induced by 
the quark contact interactions, and discuss its impact to 
the existing experimental limits set by both the ATLAS 
and CMS Collaborations. 

To compare with the experimental analyses, we con- 



sider only the quark contact interactions that are the 
products of left-handed electroweak isoscalar quark cur- 
rents which are assumed to be flavor-symmetric to avoid 
large flavor-changing neutral-current interactions [Ij. 
The effective Lagrangian can be written as 
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where ci, C2 are the Wilson coefficients, and Oi, O2 are 
the color-singlet and color-octet operators given by 
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in which c, d are the generation indices and i, j, k, I, a 
are the color indices, and T° are the Gell-Mann matrices. 
The Wilson coefficients at the scale A are conventionally 
normalized to be ci(A) = 47rcos0, C2(A) = 47rsin6', with 
< < 27r. Notice that the above operators can also 
arise from the exchange of new heavy resonances in vari- 
ous new physics models, such as Z' models Q and extra 
dimensions models 0]. Thus, our analyses are rather 
model independent and A can be identified as the effec- 
tive new physics (NP) scale. When using the above op- 
erators to calculate an observable at a scale much lower 
than A, we have to consider the QCD running effects of 
the Wilson coefficients [sl, which can be easily derived 
by solving the rcnormalization group equation using the 
one-loop anomalous dimension matrix of the operators, 
as following: 
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with r = as{fJ,R)/as{A), and 



b, = ci(A) + ^c2(A), 62 = ci(A)-^c2(A), (4) 

where iV = 3 and Cf = 4/3 for QCD, ^0 = (UN - 
2n/)/3, and /i_R is the rcnormalization scale, 71/ = 5 
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is the number of active quark flavors. In the studies 
of the ATLAS and CMS CoUaborations, they focus on 
the color-singlet operator with destructive interference, 
which corresponds to = in our analyses. Below, we 
will first discuss this case, and then extend to more gen- 
eral cases with arbitrary 9 values. 

At the LO, there are several subprocesses which con- 
tribute to the dijet production at hadron colliders in- 
duced by the NP operators we considered, including 

m'{q) qq'iq), qq' qq', qq qq{q'q'), (5) 

where q, q' could be all the quarks except the top quark. 
The NP contributions included in our calculation consist 
of two parts, the NP squared terms and the interference 
terms between the NP and the SM QCD interactions, 
which yield different kinematic distributions. We carried 
out the NLO calculations in the Feynman-'t Hooft gauge 
with dimensional regularization (DR) scheme (with naive 
75 prescription) Q in n = 4 — 2e dimensions to regularize 
all the divergences. Below, we only show the analytical 
results for the subprocess q{pi)q' {P2) — >■ q{P'i)q' (pi), since 
the similar results for other subprocesses can be obtained 
by crossing symmetry. 

First, we define the following abbreviations for the 
color structures and the matrix clement. 
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where ii_4 are the color indices of the external quarks. 
The LO scattering amplitudes induced by the NP and 
the SM QCD interactions can be separately written as 



lXo(ciCi+C2C2)/A2, 

iMQ{AiTasC2)/t, 



(7) 



where s, u are the Mandelstam variables, and we only 
keep the left-handed current product of the SM QCD am- 
plitudes here since others have no interference with the 
NP interactions. After adding the 1-loop amplitudes and 
the corresponding counterterms, we have the ultraviolet 
finite virtual amplitudes as follows. 
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We have checked that the virtual correction for the SM 
QCD contributions given in Eq. ^ agrees with the ones 
shown in Ref. @. The infrared divergences in virtual 
corrections should cancel with those in real corrections. 
As for the real corrections, we apply both the two cutoff 
method [l3| and the dipole subtraction method [ll| in our 
calculations for a cross-check. Furthermore, we do not 
include loop corrections which would induce other new 
high dimension operators and involve more uncertainties, 
and will discuss them elsewhere. 

To compare with the dijet measurements of the AT- 
LAS and CMS, we adopt the anti-fcy .iet algorithm [l^ 
with the Et recombination scheme [13| at parton level. 



We use the same parton distribution function sets [14 1 
and set both the rcnormalization and factorization scales 
to the average jet transverse momentum pT-, as in the ex- 
perimental analyses. Since our calculations are based on 
the effective field theory approach, we only consider the 
dijet invariant mas mjj up to the NP scale A. The SM 
QCD contributions are also calculated up to NLO, us- 
ing both the modified EKS code [H] and the NLOjet++ 
program [l6| for a cross-check. 

The ATLAS measurement — Following the ATLAS 
analysis, we impose the following kinematic cuts for se- 
lecting the dijet events. 

|?/jet| < 2.8, Pti{2) > 60(30)GeV, m,, > 1.2TcV, 

\yb\ - \yi + 2/21/2 < 0.75, \y*\ = \yi - y2\/2 < 1.7. (10) 

The jet radius parameter AR is chosen to be 0.6. We 
note that at parton level, the constraints on y* and rrijj 
have already set a threshold of about 150 GeV for pT2- 
To quantify the NP effects on dijet angular distribution, 
we divide the region of |y*|, between and 3.3/2, into 
11 bins with equal bin width, and define F as the ratio 
of the number of events in the first 4 bins (with cross 
section a{4:th)) to that in the total 11 bins (with cross 
section cr (tot)) [1]. InFig.^a) we show the effect of new 
physics contribution, defined through the discriminator 
A = {Fnp+sm — Fsm)/Fsm, i.e., 

^ _ Fnp+sm ^ 
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as a function of the compositeness scale A in three cases: 
(1) including only the LO results as done by CMS; (2) 
using the scaled NLO QCD results as done by ATLAS; 
(3) including the exact NLO QCD corrections presented 
in this Letter. While the scaled NLO QCD calculation 
predicts a higher value of A than the LO calculation, the 
exact NLO QCD calculation yields a much smaller value. 
To further quantify the differences, we plot in Fig. [2] the 
ratios (K-factors) of the scaled NLO (labeled as NLOl) 
and the exact NLO (labeled as NL02) predictions for A 
to the LO results as functions of A. The difference in 
the two NLO K-factors is also shown, which is generally 
larger than 30% for A larger than about 3 TeV. Using 
the scaled NLO prediction, ATLAS found the 95% CL 
exclusion limit on the quark compositeness scale is 3.4 
TeV, which corresponds to A value about 0.38. If using 
the exact NLO QCD calculation, the exclusion limit will 
reduce by about 10% to 3.1 TeV. 
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FIG. 1: A as functions of the scale A. The NLOl (2) curves 
represent the scaled (exact) NLO results. 

Fig. [Ijb) shows similar results as Fig. [T](a), but re- 
quiring the dijet invariant mass to be larger than 2 TeV 
instead of 1.2 TeV. This measurement will become feasi- 
ble when more integrated luminosity is collected at the 
LHC. Again, the scaled NLO calculation would overesti- 
mate the exclusion limit as compared to the exact NLO 
QCD calculation. For A ~ 0.38, the difference is more 
than 10%. 
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FIG. 2: The scaled NLO and the exact NLO K-factors, and 
their differences, as functions of A. 

The CMS measurement — Following the CMS analy- 



ses, we require 

hjetl < 1.3, PTjet > 50GcV, AR = 0.7, (12) 

and define the centrality ratio as the number of 
events in the inner pseudo-rapidity region I771.2I < 0.7 
(with cross section tTi„) to the one in the outer region 
0.7 < \r]i^2\ < 1-3 (with cross section aout) to quantify the 
NP effects on the dijet angular distribution, and evaluate 
R^ in each invariant mass bin from several hundred GeV 
to 3 TeV [Hj. In our calculation we use the same cuts 
and strategy as the CMS except that we only consider 
the last five invariant mass bins with equal width rang- 
ing from 1.5 TeV to 3 TeV, since the NP effects are more 
significant in this region. To quantitatively compare var- 
ious theory predictions, we introduce the discriminator 
p as the weighted average deviation of R^ from the SM 
predictions, 

5 5 
1=1 1=1 

= \l C/ {(Jin,NP+SM{i)^^ + CFout,NP+SM{i)^^), (13) 

where LOi is the inverse of the relative statistical error 
of Rjfp^sM^^)^ ^^"^ ^ integrated luminosity. In 

Fig. [IKa), we plot p as functions of the compositeness 
scale A, as predicted by the LO and the exact NLO cal- 
culations. The exact NLO QCD corrections generally re- 
duce the value of p as compared to the LO prediction. For 
completeness, the exact NLO K-factors of p arc shown in 
Fig.[3Jb). Consequently, the exact NLO QCD corrections 
modify the CMS 95% CL exclusion hmit of A from 4.0 
TeV to 3.7 TeV. 
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FIG. 3: (a) p as functions of the compositeness scale A, as 
predicted by the LO and the exact NLO calculations, (b) The 
exact NLO K-factors of p, as functions of A. 

General cases with arbitrary 9 values — It is inter- 
esting to also study the general cases in which both the 
ci(A) and C2(A) coefficients in the effective Lagrangian 
are nonzero. Here, we consider the discovery potential of 
the LHC with y/s ~ 14 TeV, and raise the dijet invariant 
mass cut to 3 TeV. In Fig. 2] we plot the discriminator A 
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as functions of 9 for different values of A together with the 
3ct error bands from the SM predictions, which arc cal- 
culated by assuming an integrated luminosity of 1 fb~^. 
We find that A varies largely from positive to negative 
values with the change of 9 for a given A value, and the 
difference between the scaled NLO and the exact NLO 
predictions can be larger than 30%. Another interesting 
observation is that there are regions of theory parameter 
space in which the NP effects on the angular distribution 
vanish, due to the cancellation of the contributions from 
the Oi and O2 operators. In that case, it will be difficult 
to study the quark compositeness scale through the dijet 
angular distribution measurement at the LHC. Finally, 
the quark contact interactions can also modify the top 
quark pair production at the Tevatron and the LHC. It 
may be possible to distinguish the effects induced from 
the two operators by studying the correlations between 
the dijet and the top quark pair production processes, 
which will be presented elsewhere. 
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FIG. 4: A as functions of 9 for different values of A. 

In conclusion, we have calculated the exact NLO QCD 
corrections to the dijet production at the LHC, induced 
by the quark contact interactions which may arise from 
the quark compositeness models or other new physics 
models. From our results, the current exclusion limits 
of the quark compositeness scale set by the ATLAS and 
CMS Collaborations shaU be lowered from 3.4 TeV to 3.1 
TeV, and 4.0 TeV to 3.7 TeV, respectively. Moreover, 
we discussed the general cases with color-octet operator 
included in the quark contact interactions and found that 
in some regions of the theory parameter space, the quark 
compositeness may become undetectable at the LHC via 



the dijet angular distribution measurement. 
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